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^S) , Abstract 

5—( ' In this article, we present a generalization of the Benamou-Brenier formula, that links the Wasserstein 



distance to an action functional related to the transport equation. On one side, we introduce a gener- 



■^^ ' alization of the Wasserstein distance for mass-varying measures. On the other side, we define an action 



functional related to the transport equation with source. 
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1 Introduction 



The problem of optimal transportation, also called Monge-Kantorovich problem, has been intensively studied 
'^k , in the mathematical community. Related to such problem, the definition of the Wasserstein distance in the 

space of probability measure has revealed itself to be a powerful tool, in particular for dealing with dynamics 
^D , of measures (like the transport PDE. see e.g. [1]). For a complete introduction to Wasserstein distances, see 

">; [6,7]. 

^^ , The main limit of this approach, at least for its application to dynamics of measures, is that the Wasser- 

^3 ' stein distance Wp(/i, z/) is defined only if the two measures /i, i/ have the same mass. For this reason, in 

^^ [5] we defined a generalized Wasserstein distance Wp''' (/i, i^), combining the standard Wasserstein and L^ 

distances. In rough words, for WS:'' {fi, v) an infinitesimal mass 6[i of [i (or v) can either be removed at cost 
a|(5/i|, or moved from /i to j/ at cost 5Wp((5/i, Ji'). This distance is a generalization of the "fiat distance", 
k^ ' One interesting field of application of this new distance is the study of transport PDEs with source, i.e. 

$_( . the following dynamics of measures: 



dt\it + V • {vt\it) = ht (1) 

Several authors have studied (1) without source term, i.e. ft, = 0, showing that it is very convenient to use 
the standard Wasserstein distance in this framework. In particular, Benamou and Brenier showed in [2] that 
there is a natural equivalence between the minimization of the action functional A [/i, v\ :— L dt (Jmd df^tl^tl'^) 
and the computation of the Wasserstein distance. Their fundamental result is recalled in Section 3. 

The main limit of the approach based on the standard Wasserstein distance is that, as explained above, 
it cannot encompass the case of a source h. Indeed, in this case the mass of the measure /it varies in time, 
hence in general Wp{^t, fJ-s) is not defined for i 7^ s. In this article, to deal with the source in (1), we focus on 
the space of Borel measures with finite mass on 'R'^ (denoted with A4), that we endow with the generalized 



'Department of Mathematical Sciences, Rutgers University - Camden, Camden, NJ. piccoliacamden. rutgers.edu 
tAix-Marseille Univ, LSIS, 13013, Marseille, France, francesco.rossiaisis.org 



Wasserstein distance WJ^' . We also denote with M.q'^ the subspace of Al of measures that are absolutely 
continuous with respect to the Lebesgue measure and with bounded support. 

Our goal here is to generalize the Benamou-Brenier formula to this setting. On one side, we define the 
generalized Wasserstein distance^, mixing creation/removal of mass and transport of mass. We first define 

7^2"' (a*, ^)= . inf, ^ a^\^Ji-lL\ + \v-u\f + h^Wi{fi,u), 
from which we have the following definition of the generalized Wasserstein distance: 

On the other side, we define a generalization of the functional A^ taking into account both the transport and 
the creation/removal of mass in (1). More precisely, we define 

B''''>[fi,v,h]:=a^( f dt( f d\ht\]] +b^ f dt( f d^lt\vt 
\Jo XJr-^ // Jo XJr'' 

Given the generalizations both for the distance and the functional, we will then prove the generalized 
Benamou-Brenier formula under the regularity hypotheses recalled in Definition 8: 

-,a.ti ^ \ ■ r I K!a,b r Ml /i is a solution of (1) with vector field v, source h 



The proof is based on two steps. We first prove (2) under some regularity assumptions for u, h, and 
approximating solutions of (1) with a sample-and-hold method. We then prove (2) in a more general setting, 
via a regularization argument. 

The structure of the paper is the following. In Section 2 we define the generalized Wasserstein distance 
and prove some useful properties. In Section 3 we recall the standard Benamou-Brenier formula and state 
the main result of the paper, that is the generalized Benamou-Brenier formula. Section 4 is devoted to the 
study of solutions of (1) under strong regularity hypotheses for u, h and to the proof of (2) in this setting. 
Finally, Section 5 contains the proof of the generalized Benamou-Brenier in the general setting, that is based 
on a regularization argument. 

2 Generalized Wasserstein distance 

In this section we define the generalized Wasserstein distance and prove some useful properties. Several 
other results, like the carachterization of the topology induced by such distance and the comparison with 
other distances, are presented in [5]. 

2.1 Notation and standard Wasserstein distance 

We first fix the notation that we use throughout the paper, and recall definitions and properties related to 
measure theory and the Wasserstein distance, like push-forward of measures 7^/i and transference plans. 

Let /i be a positive Borel measure with locally finite mass. If /ii is absolutely continuous with respect to 
fi, we write /ii <ti IJ-- If Mi(^) ^ /^(^) for a-H Borel sets, we write jii < 11. Given a measure with finite mass, 
we denote with \^\ := /i(R ) its norm. More in general, if /i = ^^ — ji^ is a signed Borel measure, we define 
\^\ :~ l/i^l + |/i^|. Such norm defines a distance in Al, that is |/i — v\. It is useful to recall that, if /ii <C /i 
and dill = f d^i with / e L^{dii), then |/xi| = / |/| d/i. 



^Observe that the definition in [5| was VFj ' {/J., v) = inf^ ,;;g^ a|/i — /i| +a\u — v\ +hWp{iJL, v). Clearly, the two definitions are 
extremely similar, and satisfy similar properties: one can indeed observe that , given the vector {a\^ — jl\+ a\v — i'\, bW2{fi-, i')) G 
M^ , the definition in [5] is the 1-norni of such vector, while the definition given in the present article is its 2-norni. 



Given two measures /i, v^ one can always write in a unique way /i = jiac + A*s such that ^ac ^ v and 
y^s -L V, i.e. there exists B such that ^s{B) = and i^(M" \ B) = 0. This is the Lebesgue's decomposition 
Theorem. Then, it exists a unique / G L^{dv) such that d/iacix) = f{x)dv{x). Such function is caUed the 
Radon-Nikodym derivative of /i with respect to v. We denote it with D^ji. For more details, see e.g. [3]. 

Given a Borel map 7 : W^ — t- M'', one can consider the following action on a measure n E Ai: 

An evident property is that the mass of /i, i.e. ^{M.'^) is identical to the mass of 7#/i. 

Given two measures fx, v with the same mass, it is thus possible to ask if there exists a 7 such that 
V = 7#/i. We say that 7 sends /i to i>. Moreover, one can add a cost to such 7, given by / [7] :— 
|/z|~^ Jjjjd |x — 7(x)|^ d^{x). This means that each infinitesimal mass (5/i is sent to Siy and that its infinitesimal 
cost is related to the p-th power of the distance between them. The, one can consider the map 7 minimizing 
such cost, if it exists. This is known as the Monge problem, stated by Monge in 1791. 

In general, this procedure works only for special /i, i> and p. Indeed, there exist simple examples of /i, v 
for which a 7 that sends fi to v does not exist. For example fi ~ 25i, v = 60 + 62 on the real line have the 
same mass, but there exists no 7 with v = 'yi^p, since 7 cannot separate masses. Moreover, one can have 
a sequence 7„ of maps such that / [7„] is a minimizing sequence, but the limit is not a map 7*. A simple 
condition that ensures the existence of a minimizing 7 is that fi and ly are absolutely continuous with respect 
to the Lebesgue measure. 

For such reason, one can generalize the problem to the following setting. Given a probability measure tt 
on R"^ X R'', one can interpret it as a method to transfer a measure p on M'' to another measure on W^ as 
follows: each infinitesimal mass on a location x is sent to a location y with a probability given by T:{x,y). 
Formally, p is sent to v if the following properties hold: 



dTT{x,-) =dp{x), \iy\ / d7r{-,y) = diy{y). (3) 

Such TT is called a transference plan from pto i^. We denote the set of such transference plans as n(/i, i^) . Since 
one usually deals with probability measures p,i', the terms \p\, |t^| are usually neglected in the literature. 
A condition equivalent to (3) is that, for all f,g E C^iMf^) it holds \p\ J^d-^^d{f{x) + g(y)) d'K{x,y) = 
Lrf f{x) dp(x) + L^ g(y) dv{y). Then, one can define a cost for tt as follows J [tt] := /ndxRd N ^ vY" dTT{x, y) 
and look for a minimizer of J in n(/i, v). Such problem is called the Monge-Kantorovich problem. 

It is important to observe that such problem is a generalization of the Monge problem. Indeed, given a 
7 sending p to v, one can define a transference plan tt = (Id x j)^p, i.e. di:{x, y) = p{W^)^^ dp{x)5y^^(x)- 
It also holds J [Id x 7] = / [7]. The main advantages of such approach are the following: first, the existence 
of at least one tt satisfying (3) is easy to check, since one can choose ti{A x B) = \p\^^ p{A)v{B)^ i.e. the 
mass from p is proportionally split to v. Moreover, a minimizer of J in n(/i, v) always exists. 

A natural space in which J is finite is the space of Borel measures with finite p-moment, that is 

MP:=\pEM\ f \x\P dp{x) < 00 



One can thus define on Ai^ the following operator between measures of the same mass, called the Wasser- 
stein distance: 

Wp{p,i^)^{\p\ min J[7r])i/P. 

7rGn(/i,i/) 

It is indeed a distance on the subspace of measures in Ai^ with a given mass, see [7]. It is easy to prove 
that Wp{kp,ki') = k^/PWp{p,v) for fc > 0, by observing that I{{kp^kv) = I{{p,v) and that J [tt] does not 
depend on the mass. We will recall some other properties all along the paper, when useful. 



2.2 Definition of the generalized Wasserstein distance 

In this section, we recall the definition of the generalized Wasserstein distance and prove some useful prop- 
erties. 

Definition 1 Let ijl,v ^ M he two measures. We define the functionals 

T^"" (m, i^) ■■= inf , , , a^ (Ia* - mI + I^^ - ^1)' + b^W^fi, u), (4) 

p,,i>^A4, \fi\—\i>\ 



W:^-\l,,v) -.^ ^T^%^). (5) 

In the following, we will also use the C^-norni for time- varying measures /x, i^ G C([0, 1], Al), that is 

d{y.,v) -.^ sup W2' {tJ't,i't)- (6) 

te[o,i] 

We now recall some properties of W2 and T"^' . We give sketch of the proofs. Detailed proofs can be 
adapted from proofs given in [5]. 

Proposition 2 The following properties hold: 

1. The infimum in (4) coincides with 

inf ami^-fi\ + \i^-i)\f + b^Wi{fL,C'), 

where we have added the constraint /i < /i, i^ < t^. 

2. The infimum in (4) is attained by some p., v. 

3. The Junctional W2' is a distance on Ai. 

4. It holds W^-" (a*, 0) < a\^i\ 

Proof. Property 1. Take a pair fl,i' E A^g^ with |/i| = |i>|, and n the transference plan realizing W2(p, v). 
Assume that fi ^ 11. Define /i' := min{l, D^/i}// and observe that, by construction fj! < fj,. Also define i'' 
the measure such that tt is a transference plan in n(/i', i^'). Then, one can prove that 

l/i - ^i'\ + \iy- u'\ < Im - Ml + 1^ - ^1, Wi{^l', v') < Wiifi, v). 

If v' < V, then the result is proven. Otherwise, applying to v' the techniques described above for jl, one 
can define v" < v and a corresponding /i" < jjl' < pL such that |/i — /i"| + |i^ — i^"| < |/i — /i| + |z/ — i>| and 
T4^|(m", v") < W^iji, i>). The resuh easily follows. 

Property 2. First use Property 1 to add the constraint /i < /i, i^ < i^. Take a sequence jl",i)" such 
that a^ (1^ - /i"| + \iy - i>"|)^ + b'^WH/l'', i>") -^ T^''' {fi, v). Since /i" < m, then the sequence has uniformly 
bounded mass, bounded by \^\. Then there exists a subsequence weakly converging to some /i*, i.e. /i" — ^ M*, 
due to weak compactness (see [3]). The same holds for v^ -^ v* . One can then prove that 

a^ (Im - m"I + I^^ - i>"l)' + &'vk|(m", ^") ^ a' (Im - m1 + I^ - ^*\f + b^wi{^l* y). 

Property 3. Symmetry is evident. 

We now prove that W2' {p-,i^) ~ implies 11 = v. Let /i, i> realize the infimum in (4). We clearly have 
l/i — /i| = |i^ — i>| = that implies ^i = jj. and v = v. Since 'W2 is a distance, then T4^2(Mi v) =[) implies /i = i>, 
that gives the result. 



We now prove the triangular inequality W2' (Mj'?) ^ ^2' if^T^) + ^2' (^tV)- Let {fl^vi) realiz- 
ing W2' (Mj '^) a-iid {v2,fl) realizing VF^' (j^, jy). Also denote with 7r^,7r^ the transference plans realizing 
W2{fl,iy^),W2{iy'^,fj), respectively. Define now D :— uim{Di^2i)^ ,lj 1)2. Define ji^fi the marginals of v with 
respect to 7r^,7r^ respectively, i.e. tt^ G n(/2, i?) and tt^ e n(i?, ^). Thus we have tt^ o tt^ e n(/2, 77). One can 
prove that 

\^i~fi\ + \ri-f]\^\^i-il\ + \il~li\^\-q-f]\ + \i)-f]\<\^i-ji\ + \v-V2\ + \il-f]\ + \v-vi\ 

W2{fi,fi) <w2{n,v) + W2{i>,n) <w2{ji,vi) + W2{y2,f])- 

Rearranging terms, one has 

T^"^^ (Ai, 77) < a" {\^l - Ml + 1^^ - ^i\f + b^Wiifi, h) + a" {\u - z^al + h - v\f + 5'T4^|(^2, ry) + 

+2a2 (1^ -li\ + \v- Vi\) {\v -V2\ + \Tl-fl\) + 252T^2(A, Vi)W2{v2,fl). 

Apply at the last two terms of the right hand side the estimate xy + zw < \/ x^ + z^^j y^ + ufl for positive 
numbers, that gives 



^2 



T'' (m, ^) < T2' (a*, y) + T2' {y. v) + 2^lT^%Iy)^T^^^) = {w^^' in, u) + W^^' (z., ,?)) ' , 

from which the triangular inequality easily follows. 

Property 4. Take the choice /i = ;/ = to estimate the infimum in T^' . D 

2.3 Topology of the generalized Wasserstein distance 

In this section we recall some useful topological results related to the metric space M when endowed with 
the generalized Wasserstein distance. 

We first define tightness in this context. 

Definition 3 A set of measures M is tight if for each e > there exists a compact K^ such that ii{W^\Ke) < 
e for all ji G M . 

We now recall the following important result about convergence with respect to the generalized Wasser- 
stein distance, see [5, Theorem 13]. 

Theorem 4 Let {/in} &e a sequence of measures in W^ , and /in, /i G M.- Then 

W2' {jJLn,^) — > is equivalent to /in ^ /^ and {/in} is tight. 

We finally recall the result of completenes, see [5, Proposition 15]. 

Proposition 5 The space Ai endowed with the distance W2' is a complete metric space. 

2.4 Estimates of generalized Wasserstein distance under flow actions 

In this section we give useful estimates for T4^^' under flow actions. Observe that these properties are proved 
for measures /i, t^ G Mq"^ only. The main reason is that, in this setting, the standard Wasserstein distance is 
realized by a diffeomorphism, that is the solution of the Monge problem, see Section 2.1. 

Proposition 6 Let Vt be a time-varying vector field, uniformly Lipschitz with respect to the space variable, 
and 0* the flow generated by v. Let L be the Lipschitz constant of v, i.e. \vt{x) ~ Vt{y)\ < L\x — y\ for all t. 
Let /i, i^ G A^g^. We have 

. w^''{^l,4>'#^l)<bt\\v\\co\^i\ 



Proof. Property 1. Let jl < ijl,v <v realize VF^' (/i, v), and T the map realizing W2{JJ., v). Then we have 

Observe now that, since 0* is a diffeomorphism and fl < ii, then \(j)*^fx — (f>^^fl\ — |/i — /i|. Observe now that 
VF2 ((?!>* #/i,</'*#i') < e2^*T/F2(/i, i>), as proved in [4, Prop. 1]. This gives the result. 
Property 2. Choose fj, — fJ. and D = 0*#/^ to estimate T2' (/i, 0*#/i). This gives 

T^^' {^iA'i^^^) <b^w^{^,,cf,'if^i) <b^\^,\ [ \x^c^'{x)fd^,<b^\^i\ [ {\\v\\cotfd^x^b^\^if\\v\\lot\ 

D 

3 Generalized Benamou-Brenier formula 

In this section we generalize the Benamou-Brenier formula (recalled below, see [2]) to W2' ■ We recall 
that the interest of such formula is to relate the Wasserstein distance between two measures fio,fJ-i to the 
minimization of the functional J \vt\'^d^t among all solutions of the linear transport equation from ^q to /ii. 
We first recall the original Benamou-Brenier formula. Observe that we deal with probability measures 
in Mf,"". 

Theorem 7 Let /ioj/^i G ^0^ where Vq'^ := A^g'^nP is the space of probability measures that are absolutely 
continuous with respect with the Lebesgue measure and with compact support. Endow Vq'^ with the weak-* 
topology. 

LetVlpo, fii) be the set of couples measure-velocity field {ii,v) := {p-tTVt)te[a.i] suchthat^^ C([0, 1] ,7'q'^), 
V e L^{dfj,tdt), Uj£[o^i]Supp(/if) is bounded, and such that they satisfy the following boundary value problem 

Define the action functional A [/i, v] := L dt [Jmd djit \vt\) on T^(/iOj A*i)- Then, it holds 

wi{^lo,^ll)^\ni{A[^l,v] |(^,w) e f(mo,mi)}- (7) 

Such result has been proven to hold also in the larger space of probability measures with finite second order 
moments, see [1]. It is also easy to prove that (7) holds for /ioj/^i G M.q'^ with the same mass to. Indeed, 
it is sufficient to use (7) for m~^ jiq., m^^fii and to observe that we have the same degree of homogeneity on 
the left and right hand sides when multiplying by a constant. 

We now prove that a similar result holds for VFj' and the transport equation with source. We first define 
the space and the functional that we study. 

Definition 8 Consider fiQ, iii G Alg^. Let V{fio,iii) be the set of triples (measure, velocity field,source term) 
{pL,v,h) := {iit,Vt,ht)te[o,i\ such that 

• /i G C([0, 1] ,M'^'^), with A^q"^ endowed with the weak-* topology; 

• V <E L"^ {d^tdt) ; 

• /i e -^"^([0, 1] , A^o^) in the sense that J^ dt (/jj^ d\ht\) < oo; 

• Ujg[o_i]Supp(/Xt) is bounded; 



they satisfy the following boundary value problem 



dtt^t + V • (vtfit) == ht, 
M|t=o = Mo, M|t=i =Mi- 



(8) 



We define the action functional 



i3"'''[/x,f,/i] ■=a^( f 



dt{ I d\ht\] ] +b' I dt{ I dfit\vt\' 



Remark 9 Observe that the conditions given above also imply that Utg[o.i]Supp(/it) C Utg[o.i]Supp(/it), and 
in particular ht have uniformly bounded support. Indeed, by contradiction, assume that Ut^[o^i]Supp{ht) (t 
UtG[o,i]Supp(/it). Looking at ft- as a functional on C^ functions, this means that there exists a function tp E 

C^ ([0, 1] X ]R'',]R) with supp(i/>) C [0, 1] x (M^ \ (Ute[o,i]Supp(^t))) and such that J^ dt J^^ dhtTpit.x) ^ 0. 

Observe now that, by construction, one has L dt J dfit{dttp + v • Vip) = 0, since ip and its derivatives are 
identically on the support of /it for each t € [0, 1]. Observe now that (/i, v, h) satisfy (1) in the weak sense. 
Choosing i/i as a test function, one has = J„ dt J„^ dhttpit^x). Contradiction. 

We now state the generalized Benamou-Brenier formula. We have 
Theorem 10 Let /io,/ii e Mq''. Then 

ini{B-^'[^i,v,h] \{^I,v,h)eV{^io,^il)}^n'''{^io,^ll)■ (9) 

It is clear the similarity between B""'' and A. In particular, the standard Benamou-Brenier formula can be 
recovered as a particular case of Theorem 10 when h = and a — > oo. 

To prove Theorem 10, we first need a series of lemmas, that we state and prove in the following section. 
We will then give the proof of Theorem 10 in Section 5. 



4 Approximation of solutions of (1) 

In this section, we give technical results that will be useful to prove Theorem 10. The idea is to approxi- 
mate solutions of (1) via an adapted sample-and-hold method, and to prove the inequality ,8"''' [/i, v, h] > 
T^' (/io,/ii) for such approximations. 

We first estimate the difference J f dji — J f dv in terms of VF^' (/i, v). 



Lemma 11 Let ^i,v e M"^" , and f e Lip(]R'',]R) n L°°(]R'', 

l,/l|oo 



fd^i 



fdv 



a.h , 



< \/2i 



) . Then 

\f\\Lip 



w^■''{^,,u). 



(10) 



Proof. Let fi < ^j-tV < v realizing W2 ' (/i, v)- We have 



f d^i 



fdv 



< 



fd{n-fl) 



fd{ji-u) 



fd{D 



< 



< ll,/l|oo|/i - Ml + \\f\\L^pW^ifl, i>) + ll/llocP - ^1, 



(11) 



where we have used that \^\ — sup {J f dfj, \ ||/||oo = 1/ and the Kantorovich- Rubinstein duality formula 
Wi(m, >^) = sup {J f d{^ — I') I ll/lliip = 1}. Recall that Wi{jl, v) < Wp{fi,i)) ioi p > 1, see e.g. [7, Sec. 
7.1.2]. Then (10) is a direct consequence of (11), by using (x + j/)^ < 2(x^ + y^). D 

We now prove the inequality B""'^ [p,, v, h] > T^' (/ig, /ii) under stronger regularity assumptions for w, h. 



Lemma 12 Let /iOjA*i G M'^f and {p,,v,h) G y(/io,A*i) such that they satisfy assumptions of Definition 8 
and additionally 

• V is uniformly L-Lipschitz with respect to x; it has C'^ -norm uniformly bounded in time, i.e. M := 
suPte[o,i] ll^tllc" < oo; 

• /i e _L°°([0, 1], Alg^), i.e. it satisfies P := s^Tp^^yQ_^ j^^ d\ht{.)\ < oo. 
Then 

B-'''[fi,v,h]>T^''>{fio,t^i). (12) 

Proof. The proof is divided into two main parts. In the first part we define ij,^", an approximation of /i 
with a variation of the sample-and-hold method, and prove appropriate convergence results of /xl*^! to /x and 
q£ ^a,6 |^^[fc] .yfc /jfcj to B"'''' [ii,v,h] for k — >■ oo. In the second part, we define a new sequence {jj}-^\v'' ,h^) 

such that S"''' jl^^\v^,h^ < S"'*" [^^''\ v'^ , h'^] and we prove that (12) holds for fi^'^^. A simple inequality 

will then give the result for general fi. 

Before the main parts of the proof, we state some simple remarks. First of all, since we deal with 
approximations of the dynamics given by v, h, then the approximated solution jv-^' could fail to be a positive 
measure for some times. Then, one needs to replace ^^^> with its positive part all along the proof. For 
simplicity of notation, this replacement is implict all along the proof. 

Second, we fix some notations that will be useful all along the proof. Given the initial datum /io, we will 
prove that all measures studied in the proof have bounded mass, and in particular |/it|, \^\ \AiA I ^ ImoI -\-P- 
We define 

m := l^ol + P- 

We also define 

r- { M L^ 

a:=V2max<^ — , - L (3 -.^ 2aP + bMm. 

Part 1: We divide this part of the proof in four parts. In the first, we define /il*^' , together with v^ ^h^ . 
In the second, we prove that ^^^> is a Cauchy sequence in a complete space, then it admits a limit /i*. In the 
third, we prove that /i* = /i, by proving that /i* is the unique solution of (1). In the fourth, we prove that 
j^a,b \^^[k] ^ yk^ ^fcj converges to B"^-^ [^l, v, /i] for fc ^ oo. 

Part 1.1: We first define /il*^!, together with w*^, h'' . The main advantage of the definion given below is 
that v'' and h'^ will never act at the same time, i.e. w^ 7^ implies h^ — Q and viceversa. This helps both in 
the definition of /i'*^! and in the computation of ,B°'^. A scheme of the evolution of the mass \ji\ \ is given in 
Figure 1. 

Fix fc e N and define At := 2^*^. We discretize the time interval [0, 1] in small intervals [nAt, (n + 1)A/:]. 
The idea of the discretization is first to divide each interval [nAt, [n + l)At] in three parts: 

[nAt, nAt + At^] , [nAt + At^, (n + l)At - At^] , [(n + l)At - At^, (n + l)At] . 

On the first part we use the negative part ft,^ of h, then the velocity v, then the poitive part /i+ of h. Clearly, 
each term must be correctly rescaled, to have /il , ,■) a ^ close to /i(„+i)At- 



We define the following vector field and the source term: 

^* - for TG (At^At-At^] 

for TG (0,At2]U (At- At^At], 



k _ j At^AP''nAt+—^i-^ir-At-) iOT T E (At\ At - At% 




for re {0,At% 
for re (At^At], 




-^t 



and 



Figure 1: Evolution of |/it|, |/i[ ' | for A: = 2 



"-nAt+T 



for T e (0, At - Ar 

j-2 



^^t^'KAt+At-Hr-iAt-At^)) fo^ re {At- At^ At]. 
We now introduce the foUowing useful notations: 

• <&f^j J 1 is the diffeomorphism corresponding to the flow generated by v'' on the time interval [^1,^2]; 

• IL\ti t2] ■~ It ^ —t ^^ ^^ ^^® mass removal given by h on the time interval [ti, ^2]; 

• Hif_^ (^1 :— J ^ h^ dt is the mass creation given by h on the time interval [ti, t2]- 

Several properties of composition of these operators will be useful in the following. First, for n integer, one 
has: 

• *[nAt,(n+2)At] ~ *[(ri+l)At,(n+2)At] ° '^[nAi,(n+l) Ai] ' 



• m 



nAt,(n+2)At] 



H 



{nAt,{n+l)At] 



H 



-rrfe 



{(n+l)At,(ri+2)At]' 



and similarly for Hi^^^ („_,_ 



[nAt,(n+2)At]- 



Second, for n even, one also has, for At = 2 '^, 

^[nAt,(n+2)At] ~ ^[nAt,(n+2)At] ' 

• ^[rAt,(«+2)At] == StnAt,(n+2)At]^ ^ud similarly for i?[„^t 
Proofs are direct consequences of the definitions and of standard properties of composition of fiows. 
We now define /il*^!. First define /ij, ' := jiq- Then, for each interval (nAt, {n + l)At], define: 



(n+2)At]- 



A*nAt+T 



l^nAt JJ-[nAt,nAt+T 



nAt,nAt+T]TrH'nAt ^[nAt,nAt+T]TrlJ-[nAt,{n+l)At] 



[k] 
-^^nAt+At-At^ 



H, 



iAt,nAt+T] 



for TG (0,At2], 

for TG (At^At- At^], 

for T e (At- At^, At]. 



It is easy to prove that /i^'^l is a solution of (1) in C([0, 1], AIq^) with velocity field v'' , source h^ :— h — h , 
and initial datum /Iq = /io. It is evident that the measure has uniformly bounded mass, in particular 

l^f'l < m for aUt e [0,1]. 



It is also easy to prove the following property for /i'" 

'lii-[nAt,nAt+T 



«l^f„At,„At+rll < aAt-'rP for r G {0,At% 



W^ 



2' (mL1,/^L1+.) < <^ h^^^^^Mm + aAtP for r G {At\ At - At\ (13) 

[6AfMm + aAtP + aAi"^(r- Ai + At2)P for r e (Af - At^, At]. 

As a consequence, one has for r G [0, At] the following estimate: 

W2'' (mLI*, A^Llt+.) < At(2aP + bMm) =: /3At. (14) 



Part 1.2: We now prove that y\'^^ is a Cauchy sequence with respect to the distance d defined in (6). 
First observe that, by substitution, the following formula holds for /^(„ 1 2)At' 

,,['^1 — ffifc 41, M <f>k M-TJ*^ I Tf*^ — 

^(n+2)At ~ ^[(n+l)At,(n+2)At]tFA^(„+l)At ~ ^[(n+l)At,(n+2)At]#=il[(n+l)At,(n+2)At] + ^ [(n+l)At,(ri+2)At] — 

= *[(n+l)At,(n+2)At]# (^*^[nAi,(«+l) Ai] # ^MnAt ~ ii[nAt,(n+l) Ai] j + ^[nAt,(n+l)At] j + 

~*^[(n+l)At,(n+2)At]#ii[(n+l)At,(n+2)At] + ^[(ri+l)At,(n+2)At] = 
= *^[nAt,(n+2)At]#A*riAt ^ *[riAt,(n+2)At] #i£[nAt,(n+l)At] + *[(n+l)At,(n+2)At] #-^[nAt,(n+l)At] + 

~*^[(n+l)At,(n+2)At]#ii[(n+l)At,(n+2)At] + -^ [(ri+l)At,(n+2)At] • 

We also decompose M[„IL2-)At ^^ using properties of composition of $'', if ,H . This gives: 

t'(n+2)At ~ ^[nAt,(n+2)At]tfM„At ^[nAt,(n+2)At] T^il[nAt,(n+l) At] + 

^*^[(n+l)Ai,(n+2)At]#*^[nAi,(n+l)At]#ii[(n+l)At,(n+2)At] + ^ [nAt,{n+l)At] + -H"[(n+l)At,(n+2) At] 

We now estimate W^' (^\^l^^,^['^^2)At)'^^^^''^^^P^'^^^^^2' lA^liAt ' mLa* ) ' i-^- the value of VF^' 
at the right extreme of the interval of discretization for fc — 1 with respect to its value at the left extreme. 
We choose n even. We have: 

'''^2 l^A^(„+2)At'A*(„+2)AtJ - '^'^2 \y^[nAt,{n+2)At]fF^^nAt i ^[nAt,{n+2)At]iFf^nAt J + 

+ VF2' (^'i'[„At,(n+2)At]#ii[nAt,(n+l)At]'*^[nAt,(n+2)At]#ii[nAt,(n+l)At]j + 

+ VF2' (^$[(„+i)At,(n+2)At]#^[nAt,(n+l)At]#ii[(n+l)At,(n+2)At]i'J'[(n+l)At,(n+2)Ai]#ii[(n+l)At,(n+2)At]j + 

+ VF2 ' (^-ff[„At,(n+l)At]' *[(n+l)At,(ri+2)At]#-^[nAt,(n+l)At] j + ^2 ' (^-^[(n+1) Ai,(«+2) At] ' -^[(n+l)At,(«+2) At] j 

< £2 T4^2°^ (^M„At 'M„Atj+0 + e2 T4^2°^ (^*J'[nAt,(n+l)At]#ii[(n+l)At,(n+2)At]'ii[(n+l)At,(n+2)At] j + 

+bAtM{AtP) + < ei^^'l^a"''' (mL'I/' , A^Llt) + 6ei^^*AtM(AtP) + bAtM{AtP). (15) 

We apply the last inequality recursively. First recall that W^' I /ig , //q 1 = and that, for a suffi- 
ciently big k, it holds e^/^^^* < 1 + 3LAt and 3LAt < 1. This gives 

(l+3LAt)"/2-l ,, ,,^, ^inLAt_^ 



< 



r (^A-t^ULlt) < 6MPAt^(2 + ^LAt) '',l,Zt-l ^ ''''''^''-^^L^ ^ 
< bMP2-^ 



L ' 
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where we have used that nAt < 1. Observe that the estimate is independent of n. Applying it recursively, 
one has 

^^2 l^PnAt'A'nAt y - ^ ^ 1-2^1/2' 

Finally, take any t € [0, 1]: for each integer fc, let rife be the biggest even number such that nfe2~'^ < t. It 
clearly holds \t - n2"''| < 2"'"+^. One has 



< 2;32" 



6MP(ei^-l) .,+i,l-2-'/2 



2-(fc+i) 



1-2-1/2 



2/32" 



where we have used (14) twice for the first term and 2'+^ times for third term. Since the estimate does not 
depend on i, one has d(/xW, /^I'^+'l) < Ci2-^ with Ci := 4/3 + 



bMP(e^^-l)^-l 1 



1-2-1/2- 



Since the estimate 



dt / (c?/z*(9Jt + vt ■ Vft) + dhtft) = 



does not depend on I and (i(/i[''l , /i^'^+'l) — > for fc — > oo, we have that /j,^''' is a Cauchy sequence. Since 
C([0, 1], Al) is complete with respect to d, then there exists a limit /x* := limj;_^oo m'*'': with ^* <E M. 

Part 1.3: We now prove that /i* == /i. We prove it by proving that it is a weak solution of (1). By 
uniqueness the result will follow. We have to prove that, for any^ ft G C|^([0, 1] x R'^), it holds 

(16) 
Observe that ^^^^ is a solution of (1) with vector field v^ , and source h^ := h — h ■ Then 

j dt n (dfii^^ddt + v^ ■ ^ft) + dhUt)\ = 0. 
We then prove (16) by proving the three following limits: 



lim 

k 



dt 



(^d^i; - d^,^^^) dtft 



0, lim 

A; 



dt[ I dfi;vfVft-d^i[''^v^-Vf 



lim 

k 



1. limt 



lim 

k 



dt (f d{ht - h''t)ft) 

L dt (J^d [dfJ-t ^ dfj\ ) dtft) =0. This is a consequence of (10). Indeed, one has 



dtffi^d^il-dii^^^^dtft 



< 



dtfe■^(^^Mf')^/2I 



\dtft\\oo \\dt.ft\\Lip 



< 



2. linit 



< d(^*,^W)V2max|^^^^,^^:^ 

a b 



0. We first fix k and Af := 2 '^, and estimate 



J^ dt [J^d d^,;vt ■ Wft - d/if'«' • v/tj 

/ dt[ d^i^^t^^VnAt+t ■ V/„At+t ) - / dr ( rf^Llt+r^nAt+r ' V/„At+r ) • (17) 

Jo \Jr'^ / Jo \Jr'' 



^The index t will be useful in the following change of variable in time. 
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Using the definition of v'^^^_^^, we have that it is for t e [0, Ai^] U (At — At^, Ai] and that for 



T € (At , At — At ] it holds f„At+T = At-2At'^ """At i ^* ,, (7- At^)- Then, after the change of variable 

At-2At^ 



T ^ t := (t - At^) .^* y , we have 



/ dr( d^,^ 



At 
At - 2At2 

At 



jfe] 

' ^At+T^nAt^ 



At-2At^ 



(r-At^) • V/„At+ 



'^* At - 2At2 



At - 2At2 7o ^^ 

To go back to (17), we estimate for each t e [0, At] the following quantity^: 

dlJ-lAt+t-^nAt+t ■ V/„At+t - / ^ ^At^Aji A.f2| At-2At2 ^^«At+t • ^/„At+At^+^'-^^f'" 



< 



dfJ.nAt+t^nAt+t ' V/„At+t - d^ 



■„At+At^ + Ai^Aii,««At+t • V/„At+t 



d/i 



[fe] 



nAt+At^ + M^^ 



WnAt+t • V/„At+t — d/J, 



nAt+At 



2^A 



t-2At2 «nAt+t • V/ 



riAt+At' 



2 I At — 2Atf 



< 



^^2'"' MnAt+t>M 



[fc] 



\/2l 



M||V/t||oo i||V/t||L,p 



(18) 



jlk] 

nAt+At 



2 I At-2At2 



M||V/„At+t-Vf 



2 I At-2At£^ II CO • 



riAt+At^ + iii^iiilf 



We estimate the first term of the right hand side of (18) via 

W2''' f AinAt+t,AiIJ'j^j^^j,_^Ai^At^J < W^''' {^^*nAt+t, IJ-nL+t) + W^''' U 



a,b ( [k] [fe] 

" ^t+t'^^nAt+At^ + ^i^^t 



We estimate W^''' ( /jLAt+t. m[^^1, , ^,2 , At-2At2 ^ ) by studying three cases: 



,W 



,w 



(a) t G [0, At ]: We observe that the evolution from MnAt+t ^'^ f^nAt+At^ ^^ given by removal of 

"] 

At+At^ + Ali^^t 

1 . We then have 



mass H'^nAt+t,nAt+At^]^ ^hile the evolution from A^iAt+At^ to A^j^i ^ , . ,2 , At-2At2 ^ is given by the 
push- forward of the diflFeomorphism $^ a, oa,2 i- 

^ ^ [«At+At2,«At+At2 + M^^tJ 



TiAa,6 / [fe] [fe] 

'''^2 I t*,iAt+t'A^ 



At+At2+ ^'-2^'^ t 



< ^2'^"(mLV.mLVa. 



TTra,6 / [fe] 



[fe] 



nAt+At' ' ^nAt+At2+ ^'-^2At2 ^ 



< It- At^lAt^ip + b 



At - 2At' 
At 



t\\v''\\com = 



It- At^|At"^P + 6Mmt. 



(19) 
(b) t e (At^, At — At^]: We observe that the evolution is given by the push- forward of the diffeomor- 



phism $^ A. ,A.2 1 ■ We have 

^ nAt+t,nAt+At'+ ^*-^,^* t 



Tya,6/ [fe] [fe] 



nAt+At'^ + ^i^^ 



t- At 



< b 



< 6|2tAt-At^ 



At - 2At'' 
At 

At - 2At2 



Ik'^llc""^ < 



Mm. 



^Here we denote with || V/t \\Lip the Lipschitz constant for V/t with respect to all t, a;-variables, even if for (18) the Lipschitz 
constant in space is needed only. 
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(c) t e [At — At^ , At] : This is similar to case 1 . We have 



nAt+At^ 



+b 



At - 2Ar 



At 



t - (At - 2At2) Mm. (20) 



We estimate the second term of the right hand side of (18) via'* 



[k] 

^nAt+At^+ ^'-/f t 



< m and 



\VfnAt+t - V,f 



.At+At^ + ^i^^ 



loo < |lV/t||L.p 



t- At 



At - 2At2 

aT 



-t 



\'^ ItWhip 



2tAt - At" 



Observe that both terms of the right hand side of (18) have a symmetry property: the value in t 
coincides with the value in At — t. 



We go back to (17) and, by using (18) and the symmetry described above, we have 



/ dti j dn*^^t^^Vr,At+t ■ V/„At+t ) - / drl C^tilfit+r^nAt+r 

JO \JR<i / Jo KJr-^ 

I / dtW^-'' (^MnAt+t , f^nAt+t) + 



^ JnAt+T 



< 



< 2v2ma}! 
+2v2max 
+2v 2max 

I.At/2 



M||V/t||oo L\\Vft\\up' 



a 


b 


M\\Vft\\oo 


i||V/jL./ 


a 


b 


M||V/t||oo 


L\\yft\\L^p^ 



/ dt(|t- At2|At~^F + 6Mmt) 
Ja 

^ /•At/2 

} / dtb\2tAt-At'^ 



dtmM\\Vft\\up 



2tAt - At^ 

At/2 




l-At/2 i'At/2 

C{At - 2At^)At^/2 + C dt |2tAt - At^] +C dt |2tAt - At 

Jo Jo 



At - 2At^ 
< Cd{fi* , Ail'^1 ) At + CAtV2 + 

At/2 




(21) 



with C = 272max| ^-^IIY'"°° ' ^^^^1^^^^^, \\Vft\\up] ■ max{l,P, 26Mto, Mm}. The estimate holds for 
k > 2, for which it holds ^^^^^^i < 2. We simply estimate (21) with Cd{fi* , fi^''^)At + CAt^/2 + 
CAt^/2 + CAt^/2 + CAt^/2 < Cd{fi*,n^''^)At + 3CAt^ by using |2tAt - At^l < At^ 
Going back to our estimate, using (17) on each interval [riAt, {n + l)At], we have 



lim 

k 





At 



dtl I diilvfVJt^dii^^^v'i-VJt 



2"-l 

< lim y 

Tl = 



d-T ( / C^AiLlt+r^nAt+r ' V/„At+r 



/ dt\ I C?t^,* At+t"«At+t • V/„At+t 

Jo VJR'' / 

< lim2'=(Cd(M*,M'''')2"'' + 3C2-3fc) = \i^Cd{ii\ ^l^^'^) = 0. 



3. linifc 



jldt{j^,d{h,^h\)ft) 



= 0. We first fix k and At := 2 , and estimate the negative part 



At 



dt 



dh„At+tfnAt+t 



/ dr [ dh^^^i^^fnAt+T 

Jo \Jwi 



(22) 



Here it is sufficient to use the Lipschitz constant in the time variable. 
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By using the definition of h and tlie cliange of variable t — At ^t, we have 



At 



dT 



dh: 



nAt+ 



r/nAt^ 



= At 



-1 



\JM'' / Ja 

Going back to (22), we estimate it by 

At 



/ ^^ ( / ^/i„At+At-ir/"At+r ) = / dt ( dh^ 

Jo KJwi / Jo Vjr'' 



At+tfn^t+t^t 



dt 







dh 



riAt+t 



ifnAt+t — InAt+tAt) 



< 



.At ^^2 

/ dtP\\ft\\up{t-tAt) = P\\ft\\up{l-At)^-. 
Jo ^ 



The same estimate can be proved for the positive part. Going back to our estimate, we have 

2'°-! , ,.At / ,. \ ,.At 



hm 

k 



f dt( f d{ht - /i^)/t) < lim V I / dt( f dKAt+tfnAt+t) - f drf f 

Jo yjM'i J ^ ~n I Jo VJR'* / Jo \JR" 



'^^nAt+r/nAt+r 



< 



2-1 ,.At 



< 



lim > / 

'^ '^^ I Jo 

n=0 ■^^ 

2*°-! 

lim > / 

*= '^^ L/n 



dt 



dt 



dh 



nAt+tJ'n-At+t 



At 



dr I d^nAt+r/riAt+r 



R'i 



dh: 



fnAt+t j - / '^'^ ( / dh^nAt+rfnAt+T j 



At+t^"'^*+* 



< 



2-2fe 

< lim2-2'=P||/jL,p(l-At)^— =0 

k Z 

We have proved that ^* is a solution of (1), with ji* G C([0, 1], Al). Observe now that /x* — /i is a solution of 
(1) with initial datum 0, vector field vt and source 0. Applying standard result of existence and uniqueness 
of solutions of (1) with zero source in C([0, 1], Al), we have /i^ = /i. Since fi e C{[0, 1], Alg^), then ^* E 
C([0,l],M[5^)too. 

Part 1.4: We now prove that i?"'*" [ij}-''\ v'^ , h'^] -^ B°-''' [iJ,,v,h] for k -^ oo. This will be a consequence 
of the two following results: 



We first estimate 



/ dti (iA^„At+t|wnAt+tN - / dr 
Jo VJR'' / Jo 

/ dt ( (iA^„At+t|wnAt+tM - / dt 
Jo VJR'' / Jo 

C^MLAt+tl"«At+tP 



"MnAt+rl^riAt+rl 
dfJ-nAt+t\VnAt+t\^ 



< 



(23) 



c?M 



.At+rl^nAt+Tl 



For the first term, observe that |wnAt+tP is bounded with bound M^ and Lipschitz in space with 
Lispchitz constant 2LM. Then we can estimate the first term with 

r* dt y2max 1^, ^1 W^^' (/.t, Mp ') < 2Mad{fi, fi^'^)At. 



For the second term, we replace v^^^_^^ with its definition and apply the change of variable t ^> t = 
{r-At^)^^§^. This gives 

^At / r \ rAt 



dt 







,W 



< 2Ma 



d^J■lAt+t\^nAt+t\' 



nAt / 

[ dW2"'(MLit+t,A^!; 



dt 





[k] 

At+At 



"^M , , o At-2At2 |"nAt+t| 
: KAt+At^+ ^* ^^ t 



< 



2 I At — 2At^ 

Kt 
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Observe that both nAt + t, nAt + Af + ^^^r^t e [nAt, {n + I) At]. Then, using (14), we have 






< 2/3At. 



Going back to (23), we estimate it with 

2Ma d{ji, ^[''1 ) Ai + AMaP At^ < Kd{n, n^''^)At + KAt^ 
defining K := max {2Ma, 4Ma/3}. 
• We now estimate 



dt 



d\ht 



dt 



d\h1 



dt 



d\ht\-d\h^\ 



We observe that, by definition of /ij , one has J 



(n+l)At 
At 



dtj^.dht = I 



(n+l)Ai 
At 



dt 



d\ht\ +d\h: 



dt Lrf dh^f. . The same 



result can be proved for the positive part. Applying it on each interval [nAt, (n + l)At], we have 



dt 



= 0. 



d\h\)) -( dti d\h^\ 
We now prove that lim^ B'^''' [/x^'^l , v'', h''] — B""'^ [/i, v, h]. Using the two previous estimates, we have 



ga,/ 



^W^^fc^/jfcj -B''^''[^x,v,h] <b^J2 K(Atd{n,tx^''^) + AtA + • a^ = 6^^ (d(^,/x['^l) + 2"^^ 



Since limfe(i(^,^W) == limfe d(^*, /x^) = 0, then lim^ S'^^'' [^W, w*^, /i*^] ^ B''-'' [fi,v,h]. 

Part 2: We now define a /il'^l, together with v'', h'', that satisfies the three following properties: 

1. /iW drives ^0 to f/^\ i.e. (/i^, w'=, /i'=) G ^(^0,^^'); 



2. it holds B' 



a.h 



j^m^yk^j^k <^a,b |-^[fc]^^fc^^fcj 



3. it holds T^''' (a*o,Mi ) < 'S"''' il^^\v^,h^ 



We divide this part in three parts. In the first, we define jv-^K In the second, we prove the properties stated 
above. In the third, we prove the final result (12). 

Part 2.1: We now define /il'^1. We first recall that {p\^\v^, h^) has a particular property: for each time 

t* it exists at most one quantity w^. , h^, , h^, that is nonzero. Moreover, if v^, is nonzero, then there exists 

a whole interval [ti,ti + At — At^] containing t* such that h^ — 0,h^ — for all t £ [ti,ti + At — At^]. 

Similar properties hold for /i , /i , with interval of length At . We call such property piecewise v-h action 
property, PVHA for short. 

We now define three transformations that preserve PVHA property. The transformation induced on the 
mass is described in Figure 2. 

Transformation DOWN V : Let {ij,,v,h) be PVHA. Let t be a time such that: wt = on the interval 
[t- At'^,t + At^]; h^ = on the interval [i- At'^,t\; hf = on the interval [t, i + At^]. Then replace 
h with h defined as follows: 

{ht for t e [0,t- At2]U (t + At^l], 

h,+^t2 forte (t-At2,i], 
ft,t^At2 for t e (t,t + At^]. 
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-^t 



Figure 2: Transformations DOWN Vf^, LEFT Ci^ and RIGHT TZi^ 



Keep V. Then the solution fi of (1) with velocity v and source h satisfies 

Mt ioTte[0,t-At'^]U{t + Af,l], 

for t e (i- Af,t\, 
for t e (t,i + Ai^]. 



Mt 



pt+Af' , ,_ 



Jt-At2 



At + /f_^.2 rfs/ii 



We use the notation V^ for this transformation applied at a certain time t, i.e. 2?t(/x) '■— ft for the 
transformation of the measure and Vidi, v, h) :~ (jl, v, h). 

Transformation LEFT C : Let (/i,w,/i) be PVHA. Let i be a time such that: h^ ^ on the interval 
\i- At + 2At^,i+ At\ ht = on the interval [t - At + 2At'^,t\; vt ^ on the interval [t,t + At^]. 
Then replace v with v defined as follows: 

'vt for t e [Oj- At + 2At^] U (t + Ai^,l], 

wt := < for t e (t- At + 2Ai^,i- Ai + 3At^], 

Ut_At2 for t e (t- At + 3Ai^,i + Ai^]. 



Also replace h with h defined as follows: 
'h7 



\^[t-At+2At^,i]) 





#h 



t+At-2At^ 



for te [0,i- At + 2At'^] U (t + At^l] 
for te {i - At + 2At^ ,i - At + 3At'^], 
for te (t- At + 3At^J+ At^], 



where ^[ti.tj] is the flow generated by v. 

Keep h^. Then the solution /t of (1) with velocity ii and source h 



/^^ 



/i satisfies 



Ht := 



A^t-At+2At2 



(*f*- 



t-At+2At-'.t~Af 



[t-At+2At'- 



aV*{sI 



t+At-2Af' 



dsh^ 



At+SAt^ 



for t e [0,i- At + 2At'^] U (i + Ai^,l], 

for t e (f- Af + 2Ai^i- Ai + 3At^], (24) 

for t e (i- At + 3At'^,i+At'^]. 



We use the notation £f for this transformation applied at a certain time t , i.e. Ctlp.) := /i for the 
transformation of the measure and Ct{^, v, h) :— {[x, i), h). 
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Transformation RIGHT TZ : Let (/i, v, h) be PVHA. Let i be a time such that: hf. = on the interval 
[t- Af,i+ At - 2Ai^]; vt = on the interval [t - At'^,t\; hf ^ on the interval [t,t + At - 2Ai^]. 
Then replace v with v defined as follows: 

forte [0,t - At'^]U {t + At - 2At^ ,1], 

for t e (f- At^,t + At-3At^], 

for t e {t + At - 3At'^ J + At - 2At'^]. 

Also replace /i+ with /i+ defined as follows: 

{h+ foriG [Oj- At2]U (t + At-2At^l], 

for te (t- At^i + At-3At^], 

*[t-t+At-2At^]#/l+_At+2At^ fo^ * ^ (^'+ ^^ - 3At2,t + At - 2At2]. 

Keep h^ . Then the solution /t of (1) with velocity v and source h := h^ — h^ satisfies 



/it for t e [0,t- At-^] U (t + At-2At^,l], 

/,, := <; $[t-t+At2]#Mt-At2 ^ for t e (t - At^, t + At - SAt^], 

fi+At-3At- + $[M-+At-2At^]# {l-tlAt-^^'^'' '^^ ^?) fo"- * ^ (t + At - SAt^, t + At - 2At2]. 



We use the notation TZ^ for this transformation applied at a certain time t, i.e. TZt{p,) := jl for the 
transformation of the measure and 7^t(/i, w, h) := (/t, v, h). 

We now fix the notation for composition of trasformations of the kind 1)^, Ci,TZf:. First, define V as 
the composition V := V^^ o 2?t^_j o . . . o V^^ o 2?^^ where 1 1 < 1 2 < . . . < t „ are all times in the set 
{0, At^, 2At^, . . . , (22*= - l)At2, 1} such that Vf can be applied. We define £, 7^ similarly. Finally, we define 
TZCD as the composition TZo CoT). 

We now apply 2?, £, TZ to /i[*=l . One can easily prove that Vi can be applied to /i^'^l for t = At, 2At, . . . , (2^*=- 
l)At, while one always has L{y}-^^) = /i^'^l and TZipP'^) ~ /i[*=l. We apply TZCD iteratively to /i^'^l. One can 
observe that, after 2*= — 1 iterations, the result is a fixed point for TZCV, i.e. 7ZC'D{7ZC'D^^ ^^^(/i^'^l)) = 
7e£X>(2'=-i)(^[fc])^ ^g jjg£jjg ~[/c] _ 7^£^'(2'=-l)(^[fc]) g^^h fixed point. 

One can observe that /i^'^l is solution of (1) for a certain v'^, h^ (depending on v^ , h^) of this kind: 

forte [0,At], 
wf = for t e [0, At] U (1 - At, 1], h't ^ {^ for t e (At, 1 - At], 




forte (1- At, 1]. 



As a consequence, we have 



'^lo~ jldrih";)- forte [0, At], 



m1" = \ ^ [A*.*] #Ma1 for t e (At, 1 - At] , (25) 

Ui'Ia* + /i-At drih^rV for t e (1 - At, 1], 

where $* is the flow generated by the vector field w* . 
Part 2.2: We now prove three properties of /il'^1. 

1. /it*=l drives /ig to /i!^ , i.e. {jl^^\v^ ,h^) <eV{iio,jj\). One can prove that, if /I e F(/io,tti ), then both 
Vfdl), Ct{ti),TZt{fi) e y(/io,t*i ), since j2 is not changed at times 0,1. Since TZCV'^^ ^^' is a finite 
composition of operators VtjCtjTZt, then ji^ ' = /iOiMi ^ Mi 7 hence (/i[*l,'D*=, h'') e V^(/io, /^i )• 
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jl^'^\v'^,h^ < B"^'^ '[p}-^\v^,h^^^. We prove it by proving that, given {^i,v,h) PVHA, it 



B'''''[}Zi{fi,v,h)]<B''^'>[fi,v,h]. 



2. It holds B"'' 
holds 
B''^" [Dtifi, V, h)] - B'^'" [fi, V, h] , S"'" [Ctifi, V, h)] < B''^' [m, V, h] 



The idea of the proof for V is to observe that adding and then removing mass has the same cost 
than removing and then adding. For C instead, observe that removing and then transporting mass is 
cheaper than transporting (a bigger mass) and then removing. For 2?, observe that transporting and 
then adding is cheaper than adding and then transporting. 

Wefirst proved"''' [^,w,/i]—S"''' jl,v,h = 0, with {il,v,h) ■.= 'Dt{fi,v,h). Observe that J^_^j2 dt(j^ad\} 
li-At-' c?i (/jid d\ht\) , then J^ dt (j^^ d\ht\) = J^ dt (J^^ d\ht\) . We then have the resuh. 

We now prove B""'^ [^, v, h] — B"^'^ jl,v,h > 0, with (/t, v, h) := >Cf (/i, w, h). Denote with where $?j ^ , 
the flow generated by v. First observe that 



dt 



t-At+2At^ 



d\ht 



dt 



t-At+2At- 
i+At^ 

dt 

t 



{^[i-At+2At'^,i]) 
■ i+At^ 



#h 



t+At-2At^ 



dK] - / dt\ I dM], (26) 



t~At+2At^ 

where we have used that mass is conserved by the push-forward of measures L^ d\^\ — L^ d\(j}4l^ij\. As 
a consequence, it holds J„ dt I L^ d\ht\ j ~ L dt (J-md d\ht\) . 

Now observe that, given two positive masses /i, v satisfying ^ < v, then L^ fd^ < L^ fdiy for / positive 
function. Using the definition of /t in (24), we have At-At+3At^ ^ l^i-At+2At^7 that also implies for 

At = ^[i-At+2At^,t-At^]#fJ-i-At+3At^ < ^[t-At+2At^ ,t-At^]#^J-i-At+2At'^ = fJ-t-At^- 

Also recalling that vt — wt-At for i G [t — At + 3At ,t + At ] , we have 

dt\ \ dflt\vt\^]< dti dMt-AtH«t-Atn = / dt{ I dHt\vt\^ 



t-At+3At^ 



t-At+2At^ 



By observing that u^ = for i e (t - At + 2Ai^, t ~ At + 3At^] and w^ = for i e (t, t + At^], we also 
have 



t+At^ 



dti / dfitlvtl'') < dti d^lt\vt\^]. (27) 

lt-At+2At^ KJm'' / Jt-At+2At^ KJm'' / 

Since {^it,vt) = ifit,vt) for t e [0,< - Ai + 2At^] U (i + Ai^, 1], then estimates (26) and (27) give the 
resuh S"'^ [^i, V, h] > B"^^ [£(/x, v,h)]. 

The proof for V is equivalent to the proof for C. 

3. It holds T^' (/xo,/i!^ 1 < B""'^ jl^^\v^,h^ . Observing the explicit structure of //['^l given in (25) one 

can observe that jl\^^ < Aq = /iO; MiIa* ^ Mi = A*i ^nd that fii_^i — ^'\^t i-At]^i^Af ^^ ^ 
consequence, one has 

2 



T-.a.() / ~[k] 



inf a^ ( 1^0 - JTiol + Imi - mil ) + fo^VF2^(mo, toi) < 



< a^ ( IaI.'^I - aSI + lAf' - A^IaJ ) + h'w^ir^,^ AriA*) 



r,2TT/2/--[fc] ,~[fe] 



At 



dt 



d{K 



k\- 



dt 



l-At 



d{hty 



+ b'wi{fif„¥^^,^^_^,^#fLf,).{2%) 
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In the first inequality we have used the choice mo — AaI' "^i ^ Ai-At ^'-' estimate the infiinum. We now 
observe that the first term of the right hand side of (28) coincides with a^ (L dt ( J^^ d\h^\]] . For the 
second term, introduce the vector field ii^ :— (1 — '2^t)vfi_2At)t+At ^-'^'^ observe that the corresponding 

[0, 



flow ^'^ satisfies $?g h^MaI ~ Mi-At- Using the standard Benamou-Brenier formula (7), it holds 



^2 [^^Av'^lAt,l~At]fF^^At 



< 



dT 



rf(«'[0,.]#Mkl)l«rl 



1-At 



{l-2At)-'^dt 



At 



-r,W 



2L~,fe|2 



<*[AM]#Mki)(l-2At)1^fl 



where we have applied the change of variable t ^ f = (1 — 2Af)r + Ai. Observing that ii^ = for 

t e [0, At] U (1 - At, 1], and that l>{=^j t]#i^A\ = Mp'' "^^ have 



r,W U-^/^lz 



W^l (aS, *fAM-A*]#AS) < (1 - 2At) £ dt (^£ d^f ' l«: 

Going back to (28), we have 

T^''{^io,^^l) < a'( / dt[ I d\h 



< I dt 
'o 



2 /•! 

2 



Part 2.3: We now prove (12). For each A: it holds T2°'^ ( Mo, Mi ) < ■8"''' mW,w'=,/i'= < S"''' [m'''', w'', /i'=] . 



6^ / rft 





dMfl|^*1^ 



ga, 



rfM^l l^f P 



mW,«',^' 



Sincelimfe |T4^2°' (mo,Mi )-W^2'' (mo, Mi) I < 1™*; T4^2°' (mi , A^i ) < linife ^(m'''', m) == 0, then limfe T2"' (mo,Mi 



,[k] 



T2 (mo,Mi)- Then 



A; \ / fc L J 
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5 Proof of the generalized Benamou-Brenier formula 

In this section we prove Theorem 10. We divide the proof in two parts. In part 1, we generalize the inequality 
T2' (Mo, Ml) ^ ^"^'^ [Mj ^7 ^]- III Part 2, we prove the converse inequality. 

Part 1: We prove that T^' (mo, Mi) ^ 'S"''' [M; ^; ''■]• We observe that we have removed hypotheses with 
respect to Lemma 12. For h, we have passed from L°° to L^ regularity. On the side of w, we have passed from 
Lipschitz continuity with respect to space and uniform boundedness to vt € L^{dtd^t)- The generalization 
for the term h is provided in Part 1.1, while the subsequent generalization for the term v is provided in Part 
1.2. 

Part 1.1: We first prove that one can pass from the case of h in L°° to the case of /i in L^. The idea is 
to define m as in the proof of Lemma 12, Part 1, and to provide similar estimates. First define 



Pn 



(n+l)2- 



dt\ht 



Now repeat Part 1 of the proof of Lemma 12. First observe that all measures have bounded mass, bounded 
by m := |mo| + /g dt\ht\- Then replace (13) and (14) by 



<apS 



w: 



{^\IL[nAt,nAt+T 
° At-2At^ ^"^ + ^P 
bAtMm + 2ap^ 



for re iO,At% 

for re (At^At- At% 

for r e (At- At2,At], 



(29) 
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and 

W^'' (mLI*, MLlt+.) < '^api + bMm^t. (30) 

With the same computations, replace (15) by 



„fc 



Observe that, given a sequence a„+i < Can + bn with ag = and C > 1, one can estimate a„ < C" ^5o + 



C" ^61 + ... + fe„_i < C" ""^ X^ILi ^«- I^ '^^'" case, we have the estimate W2' (ML+2-)At'/^(ri+2)At ) — 

elLAt^{bAtM){el^'^* + l) E^=oP^ Then, using the fact that nAt < 1 and ELo?'™ < T.T=o^ Pn = Jo dt\ht\, 
one has the fohowing estimate, independent on n: 



W^' (a^La/ULI) < 36Mei^ (/^*l^*l 



dt\ht\ I 2-K 



It is then easy to prove W2' (M„2-fciM„2-* ) — 66Me4^ (/g (it|ft.t|) 2 ^. Using (30) and choosing n even 
such that n2^'^ <t<{n + 2)2^*=, estimate 



(n+l)2' 



T4^2"'MMf' ,Aif^'') < 2a{pi+pi+^) + 2bMm2-^ + (3bMei^{ dt\ht\\ 2'^ + 2a ^ pf+' + 2'+i5Mm2 

<.(n+2)2-'= <.(n+2)2-'= 

< 2a I \ht\ dt + 022-"" + 2a \ht\dt, 



-k-l ^ 



n2-'= ^n2 



withC2 := 4bMm+6bMei^ (j^ dt\ht\Y Since |ft,t| e Li([0, 1], i?), for each e it exists V'(e) := supjgjo.i-e] /t*'^'' 1^*1 

and it satisfies V'(e) ^ for e ^ 0. Then d(/i[''l, /^I'^+'l) < 4aV'(2"''+-^) + 2"'', hence ^^''^ is a Cauchy sequence 
inC([0,l],M). 

We are now left to prove that the limit ^* — limj. ^'-''' coincides with /i. We prove it by proving that it 
is a solution of (1). Perform the same estimates as in the proof of Lemma 12, Part 1.3, with the following 
modifications: For the limit 2: 



1. In estimate (19), estimate VK^ ' [P-nkt+v P-nkt+At^ ) ^ ^* Pn^ 

2. Similarly, in estimate (20), estimate W2''' (^11*+*, Mlfit+At-At^j < ^t~^Pn', 

3. This induces a change in (21), that can be estimated with Cd(^*,^[''l)At + Cp^At + 3CAt^ 

4. This in turn induces the following change in the main estimate 
nl / /• \ / 2'=-i 



lim 

k 



J dt (j ^ dfx;vt ■ Vft - dfx^^^v^ ■ Vft) I < lim [ Cd(M*, /xW) + C2-'= J2 Pn + ^^^^'^^ I ^ 

< + Clim2"'' / \ht\dt + = 0. 
k Jo 

For the limit 3: 

1. Estimate (22) with J^' dt\hnAt+t\ \\ft\\L^p{t - tAt) < At J^' dt\hnAt+t\ \\ft\\L^p. 
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2. This in turn induces the fohowing change in the main estimate 

hm I / dt{ I d{ht - h't)ft I I < hm2-«||,M|L,p / dt\hnAt+t\ = 0. 



dt( f d{ht-h^,)f\ < \im2-''\\ ftWup [ dt\hni 
KJm'' / '' Jo 



After having proved that fj,* is a solution, one has ii* — fi hy observing that existence and uniqueness is 
guaranteed for h e L^{[0, 1],Mq'') too. 

Finally, one can easily prove B'^'^ [n'' , v'' , h''] — )■ B""'^ [^,v,h] by following the estimates of Lemma 12, 
Part 1.4. The only difference is in the estimate of the second term of (23), in which one has to replace the 
use of (14) by (30). This proves that Lemma 12 holds even with h e L^{[0, l],M^''). 

Part 1.2: We now generalize Lemma 12 to u e L^ {dt d^t) ■ The proof is a generalization of the proof of 
the Benamou-Brenier formula given in [7, Theorem 8.1], Step 2. The main idea is to introduce the variable 
rat '■= PtVt, where pt is the density of ^t, and observe that pt|wtp — \rat\'^ / pt is a convex function of pt, nit- 
Then, we write B"''^ [p,, m, h] with an abuse of notation, and observe that it is convex with respect to its 
arguments. The presence of the term h makes no difference on this point with respect to [7, Theorem 8.1], 
Step 2. 

Let {p,v,h) satisfy hypotheses of Theorem 10. Denote with fl an open bounded set containing the 
support of all pt, for t G [0, 1]. As already discussed in Remark 9, we also have that ft contains the support 
of all ht, for t E [0, 1]. Define a measure p^ E Alg^ such that its density pQ is a function in C^{M.'^) and 
satisfies Pq > 1 in fl. Define now p as the solution of dtp — h, with initial datum po ■— Pq + L dt h^ . It is 
easy to observe that p^ < pt for all t € [0, 1], then pt has a density pt satisfying po > 1 in Q,, and pt has a 
bounded support. Moreover, \pt\ < \pq\ + 2 J„ dt\ht\, then pt has bounded mass. 

For each A € [0, 1], define p^ -.^ {1 — X) pt + Xpt and m^ := (1 — X)mt. Observe that for each A it holds 
dtPt + V • m^ = h. Denote with p^ the density of p^. For simplicity of notation, fix now a certain A e (0, 1) 
and denote i] :— p'^, with density x, and n := m^. Observe that the support of n is contained in fl and that 
X is uniformly bounded from below on a neighborhood of fl. 

Define now, for each £ > 0, a regularizing kernel r^{t,x) := -j^fT'''i{x/s)r2{t/e), where ri G C^{M.'^),r2 G 
Cj?°(R+) satisfy r^ > and supp(ri) e -Bi(O), supp(r2) G (0,1). Define x^ '■— X * fe, n^ :— n * r^ 
and h^ := h * r^. Denote with -ql the measure with density Xt- It is clear that dt'q'^ + V • n^ — h^ for 
t € (e, 1 — e). Observe now that, for a sufficiently small e, it holds that x is uniformly bounded from below 
on a neighborhood of the support of n'^ . For this reason, one can define Vt '.= nf /xf , that is well-defined^ 
for all t G (e, 1 — e). Moreover, w^ is a C^ function, with compact support. Then, using the results of Step 
1.1 and time reparametrization, one has 

B'^'''[r,^,n^,h-]^B--'[rj',v^,h-] > (1 - 2e)-'TS^' (vlvl-e) >T2^' (vlvl-e) ■ 
Observe now that, again by convexity of |np/x with respect to its arguments, one has 



,e|„,e|2 



t£|2 \„. ^ ^ |2 U.|2 



\nt *rs|'^ \nt\ 



xMr = ^^ = ' < ^ * r„ 

Xt *re Xt* re Xt 

and that 

( dt [ d\hl\ ^( [ dt f d\ht\) ( f dt f dxre{t,xy\ ^ f dt f d\ht\. 

Je Jr'^ \Jo Jr'' J \Jo Jr'^ J Jo Jr'' 

This implies that B"'' [77, n, h] > B"^^ [r;^ r^^ /i^] > T^''' (77^, ?7f_^) . 

We now prove that, for e ^ 0, it holds T^' (77J, ?7f_g) -^ T^' (r^o, Vi)- Since T^' {p, v) = f VF^' {p, v) j , 

then it is equivalent to prove that T4^^' [rjl^rjo) -^ and VF^' {rj'i_^,rji) -^ 0. First observe that the 
supports of the 77^ are uniformly bounded, then the sequence 77^ is tight. Observe now that 77^ ^ 77 and, due to 
continuity of both rjl, rjt with respect to time, then 77^ -^ rjQ. Recall now that T4^^' nietrizes weak convergence 



'Here it is necessary to fix the convention 0/0 := 0. 
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for tight sequences (see Theorem 4), then W2' (?7g,?7o) ^^ 0. Similarly, one can prove W2' ('7i-e,'7i) -^ 0. 
This proves B"^'' [r], n, h] > T^^ (770, m)- 

Going back to the notation /x^, the previous result reads as B"'''^ [/i'^, m'^, /i] > T2' (/iQ,/i^) for all 
A e (0, 1). Due to convexity of yS"'**, one has 

B'''''[^l^,m^,h] < {l-X)B'''''[n,m,h] + \B''^^[fi,Q,h]^{l-\)(jdtl dnt\vt\^ + ( f dt\ht\]j + 



1 \ / /•! 

2 



^A(/ dt\ht\j<U dt j ^diit\vt\'+(l dt\ht\\\=B'''''[ii,m,h]. 

We now prove that T^'^ {n^, n^) -^ T^^ {fiQ, fii) for A -^ 0. Observe that W^''' {fiQ, Mo) ^ for A ^ 0, 
which can be easily proven by 

W^'" (m^ Mo) < W^-" ((1 - X)fio + AAo, (1 - A)mo) + W^'" ((1 - A)a.o,Mo) < AIa^oI + A|mo|, 

and by recalling that \j2q\ and |/io| are finite. Similarly, both \j2^\ and |/ii| are finite, then W^' (/i^,/ii) — > 

for A ^ 0. Then W^^' (Mi^,^^) ^ W^^' (mo,Mi), hence T^^' {^,^,^^^) ^ T^^' i^^o, ^^l)■ 
Summing up, we have 

Part 2: We now prove that inf ji?"''' [/i, w, ft,] | (/i, w, ft) e V(/io, Mi)} ^ ^2*' (MOj Mi) by giving a sequence 
{fJ' ,v^ ,h^) realizing the equality at the limit. First of all, observe that there exists^ a choice jJ,o,jJ,i such 
that 

T'2°'''(Mo,Mi) = a^dMo -MoI + Imi - Mi|)^ + &^VF2^(mo, Mi), 

and with /ig < /io,Mi < ^i- Define tjj to be the optimal map realizing T4^|(/io,/ii), that exists since /io,Mi G 
Alg^ Also define (see [7]) 

i>t{x) := {l-t)x + tip{x), i;* := (V; - Id) o V't'^ /i^ — ^j:^/io, 

and recall that (/t, v*) is the choice realizing the equality in the standard Benamou-Brenier formula (7), i.e. 

Wi{ilo,P-i)^A[tl,v*]. 

The idea of the proof is to write a dynamics first driving fiQ to flo via removal of mass, then JIq to /ii via 
push-forward of measure, and finally fli to fii with creation of mass. 
Fix an integer k and define At := 2^*^. We define v'^, h^ as follows: 

^ fO forte[0,At]U(l-At,l], ^ r-At-i(Mo-Mo) forte [0,Ai], 

^* ^= 1 (1 - 2A0-1.*. .. for t e (At, 1 - At], ^* ^= » f°^ ' ^ (^^' ^ " ^^1' 

I At"^(^i -/ii) for t e [1- At, 1]. 



l-2At 



One can prove that the expression of the solution /i of (1) with vector field v and source ft satisfies 

{(1- At"^t)/io + At"4/io) forte [0,At], 

M ^t-At^ for t e (At, 1- At], 

At "^ (f - 1 + At)^i + At"^(l-t)/ti) for t e [1- At,l]. 



^The result can be proven even without assuming the existence of fio, fii, via a double hmit and a diagonahzation argument. 
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In particular, one has {^^ ^v^ ^h^) E V{fio,iii). Similarly to proof of Lemma 12, Part 2.2, Property 3, one 
can prove that 

I dt( I d\hk\) =\no-ilo\ + \^ii-lii\, W|(Ao,Ai)= [ drf [ dflr\v;A={l-2At)[ dt ( [ diit\v. 

Jo KJm'' / Jo \Jr'^ / Jo KJr-^ 

One then has 

B'^'" [M^^^^ h'^] = a^ (Imo - MoI + Ia^i - Mil)' + b^l - 2At)-'wUfio, Ai) < (1 - 2-''+Y'T^^' (mo,Mi) • 

Passing to the limit, we have the result 

■mi{B^^''[^l,v,h] I (p,v,h)eV{fio,fii)}<lmiB^^''[fi^v\h''] <T^-''(po,f^i)- 
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